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ABSTRACT 

A significant fraction of the total photospheric light in nearby galaxy clusters is 
thought to be contained within the diffuse intracluster light (ICL), which extends 
100s of kpc from cluster cores. The study of the ICL can reveal details of the evo- 
lutionary histories and processes occurring within galaxy clusters, however since it 
has a very low surface brightness it is often difficult to detect. We present here the 
first measurements of the ICL as a fraction of total cluster light at z ~ 1 using deep 
J-band (1.2 u\m) imaging from HAWK-I on the VLT. We investigate the ICL in 6 
X-ray selected galaxy clusters at 0.8^ z ^1.2 and find that the ICL below isophotes 
fij = 22 mag/arcsec 2 constitutes 1-4% of the total cluster light within a radius R^oo- 
This is broadly consistent with simulations of the ICL at a similar redshift and when 
compared to nearby observations suggests that the fraction of the total cluster light 
that is in the ICL has increased by a factor 2-4 since z~l. We also find the fraction 
of the total cluster light contained within the Brightest Cluster Galaxy (BCG) to be 
2.0-6.3% at these redshifts, which in 5 out of 6 cases is larger than the fraction of 
the ICL component, in contrast to results from nearby clusters. This suggests that 
the evolution in cluster cores involves substantial stripping activity at late times, in 
addition to the early build up of the BCG stellar mass through merging. The presence 
of significant amounts of stellar light at large radii from these BCGs may help towards 
solving the recent disagreement between the semi-analytic model predictions of BCG 
mass growth (e.g. De Lucia & Blaziot, 2007) and the observed large masses and scale 
sizes reported for BCGs at high redshift. 

Key words: galaxies: clusters: general - galaxies: clusters: intracluster medium - 
galaxies: interactions - galaxies: evolution - galaxies: elliptical and lenticular, cD 



1 INTRODUCTION 

Dominated by dark matter and with masses rising above 
10 15 Mq, galaxy clusters are ideal regions for testing our 
understanding of astrophysical processes and revealing in- 
sights into the formation and evolution of structure in the 
Universe. The ubiquitous presence of hot gas and hundreds 
of galaxies means they can be studied out to high redshift 
in optical, X-ray and now SZ surveys, and used, not only 
as probes of the hierarchical evolution process, but also to 
measure cosmological parameters - for a recent review see 
Allen et al. (2011). 

In most clusters the brightest cluster galaxy (BCG) 
dominates the photospheric stellar light in the cluster core, 
with the BCG located close to the peak of the X-ray emis- 
sion. For example, Stott et al. (2012) find that the average 
BCG offset from the cluster X-ray centroid is only O.O3R500, 
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a result consistent with Lin & Mohr (2004), who find that 
for local massive clusters the BCG lies at an average dis- 
tance of only ~ 15 kpc from the cluster centroid - a figure 
which may rise to more like 50 kpc at 2 ^ 1 (Fassbender et 
al., 2011). 

Recent work on the mass and size evolution of BCGs 
has provided a challenge for current cosmological models 
and simulations. BCGs have been found to have signifi- 
cantly larger masses (Whiley et al., 2008; Collins et al., 
2009; Stott et al., 2010) at high redshift than is predicted 
by semi-analytic models (e.g., De Lucia and Blaizot, 2007). 
These observations have demonstrated that BCGs have al- 
ready undergone the majority of their mass evolution by 
z=l, contrary to the predictions of simulations which sug- 
gest a tripling in mass since this time. It has also recently 
been observed that the scale sizes of BCGs show very little 
increase since z^l (see Stott et al., 2011), contrary to ob- 
servations of passively evolving massive galaxies in the field 
which are observed to increase in scale size by 2-5 times 
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since z~2 without showing any significant mass increase. 
Whilst simulations of massive ellipticals show that this in- 
crease could be caused by minor mergers (e.g., Naab et al., 
2009), there is no need for as many interactions in BCGs 
since z~l to explain their modest size increase. 

The notion of the intracluster light (ICL) first arose 
from observational results that showed the surface bright- 
ness profiles of BCGs were extended, in excess of the pre- 
dicted 'classic' DeVaucouleurs (r 1 ' 4 ) surface brightness ra- 
dial profile for elliptical galaxies (Matthews et al., 1967; 
Sersic, 1968; Shombert, 1988). Estimates of the fraction of 
cluster light contained in the ICL in nearby clusters range 
from 10-50%, with the upper end set by the observations in 
the core of Coma (Bernstein et al., 1995). With the grow- 
ing realisation of the ICL's potential importance there have 
been concerted efforts on both the observational and theo- 
retical front; new deep observations of ICL have been made 
in clusters now covering a range of redshifts out to z—0.8 
(Gonzalez et al., 2005, 2007; Rudick et al., 2006; Krick et 
al., 2006, 2007; Guennou et al., 2011; Toledo et al., 2011). 
In addition, recent simulations of the ICL evolution with 
cluster evolution (e.g., Conroy et al., 2007; Murante et al., 
2007; Puchwein et al., 2010; Rudick et al., 2011) open up 
the potential of new constraints on the assembly history in 
rich clusters. 

The wide range of measured ICL fractions is at least in 
part due to the numerous methods used to treat the BCG 
extended halo and the ICL as separate components (e.g. 
Gonzalez et al., 2005; Lauer et al., 2007). Although some au- 
thors fit BCGs with a double surface brightness profile (e.g. 
Gonzalez et al., 2005; Ascaso et al., 2011) and often find an 
extended outer component many times larger than the in- 
ner component, these fits are often degenerate making these 
composite analyses difficult to interpret without dynamical 
information (Dolag et al., 2010). Furthermore the different 
definitions, assumptions and fitting procedures make litera- 
ture estimates of characteristic scale-sizes of BCGs difficult 
to compare. 

In our recent paper (Stott et al., 2011) we examine the 
surface brightness profiles of BCGs at z ~ 1 using HST data 
and see tantalising hints of an extended surface brightness 
profile (beyond a DeVaucouleurs r 1 ^ 4 profile). This has moti- 
vated us to acquire deep J-band observations with HAWK-I 
on the VLT to examine the diffuse ICL component in six 
2=0.8-1.22 clusters. We present here the first search for the 
ICL in galaxy clusters at these redshifts and compare our 
results with the cosmological simulations of the ICL by Rud- 
dick et al. (2011) who use surface brightness thresholding to 
distinguish the ICL from other galaxy components, thereby 
circumventing some of the complexities resulting from the 
different parameterisations of the extended light components 
discussed above. 

The structure of this paper is as follows: in §[2] we de- 
scribe our data, observations and reduction; in § |[we de- 
scribe the methods used in our measurements; in § [4] we 
present our results and discuss possible sources of system- 
atic error; and §s [5] and [6] contain the discussion of the re- 
sults and the conclusions drawn respectively. All magnitudes 
described in this paper are in the Vega magnitude system 
unless otherwise stated. Throughout this paper we adopt a 
ACDM cosmology with H = 70 km s^Mpc" 1 , Om = 0.3, 
Q A = 0.7. 



2 DATA 

Our sample consists of 6 of the most distant, spectroscop- 
ically confirmed, X-ray luminous galaxy clusters. The clus- 
ters are in the redshift range 0.8 ^ z ^ 1.22 and their rele- 
vant properties are listed in Table [I] The clusters were im- 
aged using HAWK-I (High Acuity Wide field K-band Im- 
ager) sited at the Nasmyth focus of UT4 (Pirard et al., 
2004). The HAWK-I camera consists of four HAWAII 2RG 
2048 x 2048 pix chips and the total field of view of HAWK-I is 
7.5x7.5 arcmin. Each detector covers 217 arcsec on each side 
and the mosaic has 15 arcsec gaps between the individual 
chips. The pixel scale of the detectors is 0.106 arcsec/pix. 

Observations took place in service mode between Febru- 
ary and April 2011, the observing nights and corresponding 
seeing conditions are listed in Table[2] The observations were 
split into 1 hour blocks to allow for ease of scheduling. Dur- 
ing each 1 hour block the cluster was centered on each of 
the 4 chips in sequence and a 10 minute exposure was taken 
on each chip, to allow sufficient depth to be reached and to 
cover the gaps, resulting in a total of 40 minutes of exposure 
per block. Each of the 10 minute exposures was made up of 
15 lots of 4 x 10 sec exposures, jittered in a random pat- 
tern inside a box of size 60 x 60 arcsec (average offset ~30 
arcsec). The large 60 x 60 arcsec random jitter was chosen 
due to the large size of the BCG halo (radii typically ~ 13 
arcsec at z=l, equivalent to ~ 100 kpc). Observations were 
performed in the J-band (central wavelength 1.2 (xm) with a 
total on source integration time of between 2.1 to 4.2 hours. 
Dark and twilight flat field images were taken at the end of 
each night of observing. 

The minimum angular distance between the Moon and 
any cluster observed was set at 40 degrees separation and 
observations were performed during grey time. The obser- 
vations were scheduled in 'Band B' sky conditions (second 
highest tier). The observations were taken at an air mass < 
2, and during the observations the seeing varied between 0.5 
and 0.9 arcsec. Photometric observing conditions were not 
requested in our application and consequently only a third of 
the observations were taken in photometric conditions, with 
the rest of the data having only thin cirrus sky transparency. 
However we obtained at least one night of photometric data 
for each cluster to which all the non-photometric data was 
normalised, and a comparison of the zero points from pho- 
tometric and non-photometric nights indicated zero point 
shifts of < 0.15 mag in the measured magnitudes of the 
stars normalised to the 2MASS (Two Micron All Sky Sur- 
vey) catalogue. 

The data were reduced using the HAWK-I data reduc- 
tion pipeline EsoRex. The data underwent basic reduction 
by calibration with a dark frame, bad pixel mask and a twi- 
light flat field frame. The basic reduced data were then sky 
subtracted using the sky subtraction package in EsoRex, 
which computes a sky frame for each exposure from a run- 
ning median of several jittered exposures. This individual 
computed sky frame is then subtracted from each individual 
exposure. The sky subtracted data were then median com- 
bined to create a fully reduced individual observing block 
of all 4 chips. An example of the raw, partially reduced and 
fully reduced data along with the twilight flat field, dark 
frame and subtracted sky frame are shown in Figure [l] 

Low surface brightness components, like the ICL, can 
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sometimes be unintentionally removed as sky during data 
reduction. In order to check that the ICL was not removed 
due to the finite jitter size during the sky subtraction we 
placed a model ICL profile into the data after dark current 
removal and twilight flat fielding but before the sky sub- 
traction was performed. The model profile used was a Sersic 
profile with Sersic index n=l, corresponding to a flat, expo- 
nential disk-like profile (see Graham & Driver (2005) for a 
detailed description of the Sersic profile), with a mean sur- 
face brightness ~21 mag/arcsec 2 and a half light radius of 
300 pixels (~250 kpc at z=l). The co-ordinate of the model 
profile was placed in each exposure so that it matched the 
jitter pattern offsets, effectively having the same World co- 
ordinates as the cluster centre. The sky subtraction was run 
over these images in the same way as for the normal data 
reduction. In several trials we found that a minimum of 83% 
of the original flux in the model profile was recovered in the 
reduced data in an annulus between 100-300 pixels (85-250 
kpc) and a minimum of 84% was recovered in an annulus 
between 300 - 500 pixels (250-425 kpc). For a visual com- 
parison we show an example of the data with a model profile 
included in Figure [2] As a sanity check we performed a sim- 
ilar analysis using a stellar profile placed near the centre of 
the clusters, recovering 95-96% of the total stellar light. We 
conclude that our running-median sky subtraction technique 
recovers the bulk of a flat extended profile typical of an ICL 
covering the core of the clusters. 

The reduction of the observations for each cluster from 
the observing blocks produced between 3 and 7 images for 
each cluster. These images were then median combined using 
the IRAF function imcombine. The surface brightness levels 
reached at lcr for each cluster are shown in Table [T] 



3 MEASURING THE INTRACLUSTER LIGHT 

Previous methods employed to estimate the ICL split 
broadly into two groups. Some adopt the approach of mea- 
suring both the BCG and ICL together. This is often done 
because of the extended nature of BCGs, thus making it 
hard to distinguish between the two components. For ex- 
ample, as mentioned in the introduction, Gonzalez et al. 
(2005; 2007) fit the surface brightness profile of the BCG 
using two DeVaucouleurs profiles. A similar method was fol- 
lowed by Toledo et al. (2011), where all the galaxies in the 
cluster were masked out except for the BCG and the sur- 
face brightness profile of the extended BCG component was 
measured. This kind of measurement naturally provides es- 
timates of the fraction of cluster light which is contained in 
the BCG+ICL. 

The other group aims to determine the ICL contribution 
by removing light from galaxies and objects, including the 
BCG, above a surface brightness limit, and then integrating 
the remaining surface brightness across the whole cluster 
(Krick et al., 2006; 2007). Feldmeier et al. (2004) measure 
the flux in the ICL by masking out all the flux below a series 
of surface brightness levels and subtracting the remainder 
from the total flux in the cluster. 

Our method for measuring the ICL is of the second type 
and was motivated so as to follow as closely as possible the 
predictions of Rudick et al. (2011) (hereafter Rll), who use 
GADGET- 2 to simulate the evolution of 6 galaxy clusters with 



Table 2. Observing nights and seeing conditions for our clusters 
observed with HAWK-I. 



Cluster 


Observing 


Date 


Seeing 




Block 




f 'i r/*cf>r i 

1 ell L&UL I 


CL J0152 


1 


14/09/10 


0.87 




2 


14/09/10 


0.57 




3 


14/09/10 


0.89 




1 


io t r\r\ t "i r\ 

18/09/10 


0.56 




5 


26/10/10 


0.80 


XLSS J0223 


1 


16/10/10 


0.58 




2 


16/10/10 


0.59 




3 


16/10/10 


0.56 




4 


16/10/10 


0.65 




5 


11/12/10 


0.88 


XLSS J0224 


1 


29/10/10 


0.74 




2 


29/10/10 


0.69 




3 


11/12/10 


0.68 


RCS J0439 


1 


14/09/10 


0.70 




2 


05/10/10 


0.75 




3 


16/10/10 


0.57 




4 


28/12/10 


0.47 




5 


22/02/11 


0.73 


MS 1054 


1 


06/01/11 


0.46 




2 


22/02/11 


0.81 




3 


22/02/11 


0.57 




4 


22/02/11 


0.56 




5 


22/02/11 


0.60 




6 


05/03/11 


0.69 




7 


05/03/11 


0.62 


RDCS J1317 


1 


22/02/11 


0.82 




2 


05/03/11 


0.65 




3 


12/03/11 


0.51 




4 


12/03/11 


0.79 



simulated cluster masses ~ 10 Mq, similar to the masses 
of our clusters (Table [I]). With special emphasis on the dy- 
namics of the ICL and galactic outskirts their method was 
first presented in Rudick et al. (2006), then overviewed and 
updated in Rll. They predict the light from their simulated 
galaxy clusters below a series of surface brightness thresh- 
olds calculated between redshifts 0-2, covering the redshift 
range of our clusters. Hence we choose to measure the light 
in the galaxy clusters below a series of 7 surface brightness 
levels and count all the light below these levels as the ICL. 
Although this definition of the ICL does result in some loss 
of information regarding the contribution of the extremities 
of galaxies to the ICL, it is operationally effective requiring 
no separation of galaxies from ICL and makes comparison 
with the evolutionary predictions very straightforward. 

Before any analysis was carried out all the stars and 
point sources in the cluster images were masked out. Point 
sources were identified as having a flux FWHM consistent 
with the average seeing for the combined image of that clus- 
ter (0.5 — 0.9 arcsec). We also masked out any objects with 
integrated magnitudes greater than that of the BCG to at- 
tempt to eliminate contamination from foreground galaxies 
(a more in depth discussion on the effects of contamina- 
tion is found in § |4.1| |. We used SExtractor (Bertin & 
Arnouts, 1996) to measure the total isophotal flux in all the 
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Figure 1. An example of the steps in the data reduction process. Top row: the raw data (top left), the twilight flat field (top centre), 
the dark frame (top right). Bottom row: the data after basic dark frame removal and twilight flat fielding (bottom left), the resultant 
sky calculated from the running-median technique that is subtracted (bottom centre) and the data after all the reduction described 
above (bottom right). The cluster shown is RDCS 1317 when observed for a quarter of an observing block on Chip 1. The size of each 
individual image is 217" x 217", equivalent to 1736 x 1736 kpc at z=l. 
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Figure 2. A model surface brightness profile is added into the data after it has undergone basic calibration by dark frame and twilight 
fiat subtraction (left). The sky that is subtracted (centre) shows no visible evidence of the profile (despite an rms pixel-to-pixel variation 
of 0.045 counts compared to the flux of the model ICL ~ 1.0), while the model surface brightness profile is seen almost fully recovered 
in the sky subtracted data (right) at the position indicted. The cluster shown is MS 1054 when observed for a quarter of an observing 
block on Chip 2. The size of each individual image is 217" x 217". 
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Table 1. Clusters observed with HAWK-I, all magnitudes are in the Vega system. 



Name 


RA 


Dec 


Redshift 


Integration time 


Limiting magnitude 


Cluster mass a 










Hours 


J mag/arcsec 2 


x 1O 14 M 


CL J0152 


01 52 41 


-13 57 45 


0.83 


3.5 


24.6 ± 0.2 


4.5 ± 2.7 


XLSS J0223 


02 23 03 


-04 36 22 


1.22 


3.5 


24.5 ± 0.1 


1.8 ± 0.9 


XLSS J0224 


02 24 00 


-03 25 34 


0.81 


2.1 


24.2 ± 0.1 


2.3 ± 1.4 


RCS J0439 


04 39 38 


-29 04 55 


0.95 


3.5 


24.5 ± 0.1 


0.5 ± 0.4 


MS 1054 


10 57 00 


-03 37 27 


0.80 


4.2 


24.3 ± 0.2 


8.5 ± 4.9 


RDCS J1317 


13 17 21 


29 11 18 


0.81 


2.8 


24.2 ± 0.1 


2.7 ± 2.9 



a From Stott et al. (2010) (see references within) 

objects in the cluster (including the ICL) above the back- 
ground, out to a cluster radius of -R500, corresponding to 
the radius at which the mean cluster density is 500 times 
the critical density. ICL 'images' were made from the clus- 
ter data by placing an upper limit on the flux in the im- 
ages corresponding to surface brightnesses in the J-band of: 
fij = 18, 19, 20, 20.5, 21, 21.5, 22 mag/arcsec 2 . This was done 
using the imarith function in IRAF - the total isopho- 
tal flux in these images above the background level was 
then measured using SExtractor. The background used 
by SExtractor was measured globally (BACKPHOTO 
type GLOBAL), so the sky was taken over the whole im- 
age and local variations were not subtracted. The spatial 
scale over which SExtractor identified the background is 
set by the background mesh size parameter (BACK_SIZE). 
In our analysis this was set to be larger than the average size 
of the objects in the images at 50 pixel 2 and was smoothed 
over 5 meshes. The fraction of light contained in the ICL is 
thus defined as: "the measured flux below a given surface 
brightness limit divided by the total flux in the cluster out 
to ,R 5 oo" ■ 

The fraction of the total cluster light contained in the 
BCG was also measured. The surface brightness profiles for 
the BCGs were measured using the ellipse routine in the 
IRAF STSDAS package. To do this, first all other stellar and 
non-stellar objects in the cluster area are masked out to pre- 
vent contamination of the surface brightness profile. The flux 
is measured in ellipses at increasing radii from the centre of 
the BCG along its major axis. The ID profile produced from 
this is then fitted with a model surface brightness profile us- 
ing a chi-squared minimisation technique, similar to that de- 
scribed in Stott et al., 2011. The only significant difference 
between the method described in Stott et al. (2011) and 
the method used here is that the routine used here allows 
for the fitting of a point spread function (PSF) convolved 
profile, this allows us to fit the inner regions of the surface 
brightness profile - this correction was not needed for the 
HST data considered in Stott et al. (2011) due to the much 
smaller PSF size of HST. We fit the BCGs with DeVau- 
couleurs r 1 ^ 4 surface brightness profiles, which are shown in 
the Appendix. 

To calculate flux contained in the BCG the best sur- 
face brightness profile fit is transformed into a 2D model of 
the BCG, with position angle and ellipticity matching that 
in the original data. To maintain consistency with measure- 
ment of the ICL, the flux in the model BCG is also measured 
using SExtractor out to the same distance from the clus- 



ter centre (-R500) and the flux in the BCG is expressed as 
the fraction of the total flux in the cluster (see Table [3). 

To provide a conservative calculation of the errors the 
data were re-stacked so that a stacked image containing the 
cluster was made for each detector chip from all of the nights 
observed, i.e., all the observations of the cluster when it is 
positioned on chip 1 from each of 3, 4, 5 or 7 observing blocks 
(see Table [2} were stacked together. The error quoted was 
calculated from the variation between the measurement of 
the ICL on each of the 4 individual chips. Each individual 
chip stack reaches a shallower depth than the full stack of 
all the chips so we do not quote the ICL levels measured on 
these, merely the distribution of values between the mea- 
surements on each chip. The actual errors quoted are the 
equivalent 68% percentile values (1 standard deviation for a 
Gaussian) from a Student's t-distribution of the ICL values 
measured on the 4 chips (Student's t-distribution being a 
more robust measure of the error for small number statis- 
tics, 4 in this case). 



4 RESULTS 

Our results are summarised in Table [3] where we list the 
percentage of light in the ICL at 7 different surface bright- 
ness limits in the J window. We make significant detections 
of the ICL in these high-z clusters down to isophotal lev- 
els of fj,j = 22 mag/arcsec 2 . At this limit we estimate the 
ICL to contain ~ 1 - 4% of the total light of the cluster. 
Illustrations of the measured ICL component at some of the 
surface brightness levels measured are presented in Figures[4] 
-[9] and clearly show that the detected extended light follows 
the galaxy light very closely. We return to this point in § [5] 

Also shown in Table [3] is the percentage of the J-band 
light from the BCG estimated from the best DeVaucouleurs 
model fit, described above. For comparison with earlier work 
(e.g., Whiley et al., 2008; Aragon-Salamanca et al., 1998), 
we also measure the fraction of the J-band light in the BCG 
within a 50 kpc radius centred on the BCG in the reduced 
cluster images. 

In order to compare our results with RU we transform 
their rest frame V-band isophotes (fj,v,rest) to predict our 
observed J band (jJ,j tPr ed) using the expression; 

fij.pred(z) = ^V,obs{z) + 2.5Zoff(l + zf — Kjv{z). (1) 

Here 2.5 log(l+z) 4 is the familiar surface brightness cos- 
mological dimming term and Kjv takes into account both 
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Figure 3. An illustration of the types of measurement of flux. Above are two elliptical galaxies in the cluster CI J0152 (left). The surface 
brightness is measured above the flux level equivalent to fij = 20 mag/arcsec 2 . Centre is a pixel map of all the pixels with flux above 
this level, in our measurement of the flux in the ICL this is the method we have used. Right is an illustration of the elliptical apertures 
generated by S EXTRACTOR corresponding to this isophotal level. The size of these images is 6.4"x6.4" (50x50 kpc). 



Table 3. ICL percentages of the total cluster light in our cluster sample measured above surface brightness limits in the J-band. We also 
show percentages of the total cluster light contained in the BCG using two different measurements for comparison with previous studies. 



% cluster light at surface brightness limit /ij % cluster light in BCG 



Cluster 


18 


19 


20 


20.5 


21 


21.5 


22 


DcVaucoulcurs 


50 kpc 


















model 


aperture 


CL J0152 


84.7 ± 5.2 


73.6 ± 35.2 


51.6 ± 3.1 


32.5 ± 2.9 


13.6 ± 1.9 


5.4 ± 1.6 


2.7 ± 0.4 


3.8 ± 0.6 


2.0 ± 0.2 


XLSS J0223 


90.9 ± 4.2 


82.0 ± 21.9 


55.6 ± 2.5 


29.9 ± 1.0 


10.1 ± 0.8 


3.6 ± 0.8 


2.5 ± 0.4 


2.0 ± 0.5 


1.6 ± 0.2 


XLSS J0224 


89.2 ± 27.4 


77.6 ± 12.4 


34.7 ± 11.6 


8.7 ± 2.0 


5.3 ± 0.9 


2.6 ± 0.8 


1.3 ± 0.4 


6.3 ± 0.7 


3.7 ± 0.3 


RCS J0439 


80.0 ± 5.3 


70.1 ± 13.0 


47.2 ± 3.5 


25.7 ± 3.7 


9.3 ± 1.4 


2.8 ± 2.6 


1.5 ± 0.7 


3.3 ± 0.4 


1.3 ± 0.2 


MS 1054 


69.2 ± 7.9 


60.2 ± 15.1 


41.1 ± 4.2 


23.3 ± 1.5 


11.3 ± 2.7 


6.2 ± 0.8 


3.8 ± 0.2 


4.7 ± 0.7 


1.8 ± 0.3 


RDCS J1317 


78.5 ± 14.9 


66.0 ± 12.6 


32.4 ± 1.0 


12.4 ± 5.6 


4.6 ± 1.1 


2.5 ± 0.4 


1.5 ± 0.4 


2.5 ± 0.5 


1.3 ± 0.2 



the fc-correction for accurate waveband comparison appro- 
priate to the adopted Vega magnitude system (see Hogg 
et al., 2002) and the evolution of the stellar population. For 
the latter, Rll adopt a non-evolving fixed-age model of ~ 10 
Gyr at z = 0.8, 1.0, 1.2, which, as they admit, does not allow 
direct comparison with real clusters. However, to transform 
their predictions we also assume that the ICL is made up 
of an old stellar population, similar to the BCG's. This as- 
sumption of an old stellar population is justified based on 
the findings of other studies. Krick et al. (2006; 2007) ob- 
served the ICL in 10 clusters at 0.05 < z < 0.3 and find 
the B — r colour to be consistent with simple passive stellar 
evolution, lying on the same red sequence as the cluster el- 
liptical galaxies and find the stellar populations in the ICL 
to have metallicities between Zq and 2.5Z© based on these 
colours. Further evidence that the two stellar populations 
have the same origin comes from the lack of any strong ra- 
dial gradients in the B — V colours of M87 in Virgo (Rudick 
et al. 2010). Finally, Zibetti et al. (2005) analyse the stack of 
683 clusters between z = 0.2 - 0.3 from SDSS-DR1 finding 
that the ,r— and i— band colours of the BCG and ICL 
are identical within the statistical uncertainties. 

We transform the fjby surface brightness values from 
the simulations of Rll to J-band surface brightness using 
the Bruzual & Chariot (2003) GALAXEV stellar population 
synthesis models. We used values for the formation redshift 
(zf) of 3 and 5 and metallicity of the stellar population of 



Zq and 2.5Zq. The range of Zf between 3 and 5 is similar to 
the stellar ages found for high redshift BCGs (see Collins et 
al., 2009; Stott et al., 2010). The corrections for these two 
values of formation redshift and metallicity are listed in Ta- 
ble]^ Using the Bruzual & Chariot models over this range 
of parameters gives rise to a mass-to-light ratio (M/L) ~ 4 
- 6, which is consistent with the M/L of 5 used in the stellar 
evolution modelling of the ICL by Rll and characteristic of 
the older stellar populations expected in galaxy clusters. 

In Figure [lO] we plot the average of our results at our 
5 faintest fij levels vs the ICL fraction along with the 
Hj.pred predictions from the simulations of Rll at z = 
0.8, 1.0, 1.2. The simulations of Rll span a range of py=25- 
27 mag/arcsec 2 , corresponding to predicted near-infrared 
isophotes of (ij, pre d 24 - 26.5 mag/arcsec 2 , with only slight 
variations depending on the stellar population applied (see 
Table |4|. The ICL fractions we find (~ 1 - 4%) fall below 
the Rll predictions at 2-3 magnitudes brighter than the 
Rll surface brightness limit. 



4.1 Contamination 

One of the potential problems in measuring ICL fractions, 
particularly for clusters at such high redshift, is contami- 
nation from non-cluster galaxies. As mentioned in § [3] we 
mask out any objects brighter than the BCG in an attempt 
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Figure 4. The combined VLT image of CL J0152 (left). Proceeding images show the cluster with all light above given surface brightness 
levels blacked out to illustrate the light which is measured as ICL. The surface brightness limits correspond to fij = 19 (centre) and 21 
(right) mag/arcsec 2 . The size of the images shown here is 120"xll0" (960 x 900 kpc). 




Figure 5. As Figure[I]but for XLSS J0223. 



Table 4. Total corrections applied for surface brightness dim- 
ming, fc-correction and evolution correction for ICL at z=0.8, 1, 
1.22, redshifts matching simulations we are comparing to, using 
stellar population models of differing redshift formation (zf) and 
metallicity(Z). 



Redshift 0.8 1 1.2 



2.51og(l+z) 3 








1.91 


2.26 


2.57 


2.51og(l+z) 4 








2.55 


3.01 


3.42 


Evolution 


z f 


=3 


Z 


0.77 


0.99 


1.27 


correction 


z f 


=3 


2.5Z Q 


0.80 


1.07 


1.39 




z f 


=5 


Z 


0.68 


0.83 


1.05 




z f 


=5 


2.5Z 


0.72 


0.85 


1.14 


K (colour) 


z f 


=3 


Z 


2.0216 


1.9369 


1.8527 


correction 


Zf 


=3 


2.5Z 


2.1375 


1.9864 


1.8660 




z f 


=5 


Z 


2.0322 


1.9409 


1.8575 




z f 


=5 


2.5Z 


2.1592 


1.9994 


1.8777 


Total 


Zf 


=3 


Z 


-0.8816 


-0.6669 


-0.5527 


correction 


z f 


=3 


2.5Z 


-1.0275 


-0.7964 


-0.6860 




z f 


=5 


Z 


-0.8022 


-0.5109 


-0.3375 




z f 


=5 


2.5Z 


-0.9692 


-0.5894 


-0.4477 



to correct for contamination. Not correcting for this at all 
certainly affects our results; for example, if we measure the 
ICL fraction with just the point sources masked, leaving all 
the cluster and non-cluster galaxies in, the resulting frac- 
tions ICL at fij = 22 mag/ arcsec 2 reduce to about one half 
those given in Table [3] 

Attempts to correct for this contamination reported in 
the literature are somewhat patchy and dependent on other 
data available. Gonzalez et al. (2005) mask out point sources 
but do not discuss contamination from non-cluster members 
for their z ^ 0.13 sample. Similarly, using various samples, 
Feldmeier et al. (2004), Zibetti et al. (2005) and Krick et al. 
(2006) conclude between them that contamination by non- 
cluster galaxies is not significant out to z ~ 0.3. A more ro- 
bust approach is taken by Toledo et al. (2011), who account 
for non-member contamination by calculating the luminos- 
ity function of the cluster from spectroscopically confirmed 
members then integrating the fitted luminosity function out 
to faint magnitudes to find the total luminosity of the clus- 
ter. 

In the same vein and in the absence of extensive red- 
shift information for most of our sample, we estimate the 
effect of non-member galaxy contamination using the ex- 
tensive spectroscopic survey of galaxies belonging to one of 
our clusters, CL J0152 (Demarco et al. 2010). These authors 
provide spectroscopic confirmation of cluster membership 
for 134 galaxies brighter than K S {AB) ~ 23 (or K s = 21.9 
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Figure 7. As Figure [I] but for RCS J0439. 



in the Vega system), equivalent to J ~ 23.7 for an old stellar 
population at z = 0.8. To test for the effect of contamina- 
tion we masked out all the objects in our data within -R500 of 
the cluster centre which were not confirmed cluster members 
according to the Demarco et al. (2010) spectroscopic defi- 
nition. Performing the ICL analysis again in an otherwise 
identical way, resulted in a 20% increase in the ICL fraction 
(e.g. 3.2 ± 0.4 % c.f. 2.7 ± 0.4 % at fij = 22 mag/arcsec 2 ). 
This is a strong test of the contamination of the total clus- 
ter light measured, rather than the total ICL. Although this 
tests affects the total cluster flux more strongly, we expect 
some of the flux in the ICL to be subtracted as well when 
non-cluster galaxies are removed. 

Some increase is to be expected of course as bright clus- 
ter galaxies not targeted by Demarco et al. (2010) were also 
excluded in this test but in any case within the errors there 
is convincing evidence that contamination from non-cluster 
members is not a serious concern, as is illustrated by this 
example. 



ure 11 we show the resulting sky background values (in nor- 



malised counts) on each of the final mosaiced images for 
our six clusters, calculated using the modal value of each 
slice. Fitting a simple straight line to these data, the largest 
linear gradient on the chip in either the x or y direction oc- 
curs in the mosaiced image of XLSS J0223, which shows a 
change of 0.1 in the normalised counts over the central 1200 
pixels, corresponding approximately to the region within 
7?50o- Since the standard deviation in the background noise 
ranges from 0.26-0.34, we argue that the large-scale flat- 
field gradients in our sample are not significant and we do 
not attempt to correct for them or add a separate error com- 
ponent to account for them (for example, Feldmeier et al. 
2004). In the same units as Figure 11 the ICL is clearly 
detected at 2 — 3<r with the six measured ICL values rang- 
ing from 0.83-1.1, corresponding to surface brightnesses of 
fij = 23.48 - 22.92 mag/arcsec 2 . 



4.2 Flat Fielding 

In order to recover accurate ICL values it is necessary to 
control the flat-fielding errors, particularly on large scales, 
as these are often the dominant source of uncertainty. The 
design of the 15 jittered exposures making up each observ- 
ing block and described in § [2] enables us to use the data 
reduction pipeline EsoRex to compute the sky background 
and thus define a flat-field pattern for each cluster. In Fig- 



5 DISCUSSION 

At a surface brightness threshold limit of fij = 22 
mag/arcsec 2 we estimate that the ICL constitutes 1 - 4% of 
the total light of the clusters in our sample. These values are 
similar within the errors to those found in the simulations of 
Rll, who find ICL fractions 4 - 12%, but are measured at a 
significantly brighter isophote, with Rll reaching (Xj = 24- 
27 mag/arcsec 2 - an isophotal level at least 2 magnitudes 
fainter than our data. In Figure [To] we also see a steeper 
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Figure 8. As Figure [3] but for MS 1054. 




Figure 9. As Figure [3] but for RDCS J1317. 



Simulations of R 1 1 : 

z=1 
z=1.2 
This work: 
Observed ICL 
Observed BCG 




(li mag/ arcsec* 

Figure 10. ICL fraction vs surface brightness threshold in J at 
different redshifts. Coloured lines show the predictions of Rudick 
et al. (2011) at 2=0.8 (green), 2=1 (red), 2=1.2 (blue). The aver- 
age BCG fraction of the total cluster light at each surface bright- 
ness threshold is shown in dark blue (purple). 



relation between the ICL fraction and surface brightness 
limit compared to Rll. This may be evidence that there 
is less extended ICL than is predicted by the simulations 
at this redshift. A glance at Figures [4] - [9] indicates that 
the ICL is concentrated close to the cluster galaxies in the 
core rather than smoothly distributed. It may be that an 



ZL J0152, a = 0.29 



XLSS J0223, = 0.34 



XLSS J0224, o=0.29 



RCS J0439, 0=0.33 




200 400 600 800 1000 

pix along y axis 

Figure 11. The measured flat fielded background on our clusters 
along the y-axis of our cluster images. The mode was taken for 
each slice of image x-axis at the image y-values shown. The a 
values quoted are the standard deviation of the noise shown on 
these plots. The measurements have been offset along the y-axis 
of this plot for ease of viewing. A very similar plot is seen for slices 
along the x-axis of the images and is thus not included here. 



extended component is below our sensitivity threshold, al- 
though it is known that measured ICL components closely 
follow the galaxy distribution in a number of low redshift 
clusters (Feldmeier et al. 2004). 
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As described in § [T] due to the different definitions em- 
ployed and the range of methods used to measure the ICL 
direct comparisons with other observations are difficult. Also 
there are a wide range of ICL fractions reported between dif- 
ferent authors and within individual samples, nonetheless we 
point out some general trends here for work that uses a simi- 
lar analysis to our own. For example, Feldmeier et al. (2004) 
measure the ICL in 4 non-cD clusters in the redshift range 
2=0.15 - 0.2, using a similar surface brightness threshold 
technique, finding that the ICL below fj,y = 26 mag/arcsec 2 
contains 16 — 28% of the total cluster light, falling to 3 — 9% 
below nv = 27 mag/arcsec 2 . Krick et al. (2006, 2007) mea- 
sure the ICL in 10 clusters at z — 0.05-0.3 and find that 
the ICL below nv = 26 mag/arcsec 2 contains 4-21% of the 
total cluster light. Zibetti et al. (2005) use the SDSS to mea- 
sure the ICL in 683 clusters at 2=0.2-0.3 by plotting sur- 
face brightness profiles of clusters out to 700 kpc, defining 
the ICL to be all the light below fi r = 25 mag/arcsec 2 and 
masking out light above this level. Within 500 kpc they find 
that the ICL contains 10.9 ± 5.0% of the total cluster light, 
with the ICL being more centrally concentrated than the 
galaxies and the ICL fraction decreasing slowly with radius. 
The ICL fractions we measure in our high-z clusters are 
generally lower than those found by other studies at low-z 
where similar ICL measurement methods have been used, 
which suggests a growth in the fraction of ICL of perhaps a 
factor of 2 - 4 since z=l. 

Insight into the mass assembly process at the centres 
of clusters comes from comparing the relative contributions 
from the BCG and ICL to the total light. From Table fj the 
fraction of cluster light contributed by both the ICL and the 
BCG, measured as described in ij3] ranges from 4.0 - 8.5%. 
This is smaller than the fractions reported by previous stud- 
ies at low redshifts (e.g., Gonzalez et al., 2005, 2007; Zibetti 
et al., 2005; Toledo et al., 2011) and the predictions (Con- 
roy et al., 2007), which find that the BCG+ICL constitutes 
as much as ~ 33 - 89% of the total cluster light. However, 
whereas the majority of the BCG+ICL (up to 80%) is in 
the ICL at low redshift, at z=l we are seeing the biggest 
contribution to the BCG+ICL component coming from the 
BCG, constituting ~ 60% of the BCG+ICL compared to 
a contribution of 40% from the ICL. These results can be 
understood as a direct consequence of the stripping of sur- 
rounding galaxies in the cluster cores thus increasing the 
ICL component over time and further indicates that activity 
between galaxies in the centres of rich clusters may involve 
more stripping than merging, with matter from interactions 
ending up in the ICL or cD halo rather than centrally in the 
BCG. 

If the ICL does indeed show significant evolution over 
the same timescale where the BCG does not, it would be 
further evidence that these two components have separate 
evolutionary paths and separate origins. If the two distinct 
merging and stripping histories are correct, then one might 
well expect fossil groups, as the most dynamically evolved 
clusters, to exhibit large amounts of ICL surrounding their 
central high mass BCGs (see Harrison et al., 2012). 

The quantity of extended halo light surrounding the 
central BCGs in clusters may go some way towards explain- 
ing why previous estimates of the stellar mass in BCGs at 
high redshift using i<"-band light as a mass proxy are sig- 
nificantly larger than predicted from semi-analytic models 



(Aragon-Salamanca et al., 1998; Whiley et al., 2008; Collins 
et al., 2009; Stott et al., 2010; 2011). This idea originates 
from a comment made in Whiley et al. (2008), that despite 
large BCG stellar masses at z~ 1, much stellar mass growth 
may be taking place outside of their fixed aperture magni- 
tude measurements, made with a metric circular aperture of 
37 kpc diameter; similarly the aperture sizes in the Collins 
et al. (2009) and Stott et al. (2010) papers range from 11- 
30 kpc diameter. Testing this result we find (Table [3]l that 
the contribution from the stellar light within a 50 kpc di- 
ameter aperture centred on the BCG is similar on average 
to the ICL contribution at /ij=22 mag/arcsec 2 measured 
outside that region, and constitutes about 50% of the total 
estimated BCG light based on the 2D DeVaucouleurs model 
profile extrapolated to Rsoo- These simple comparisons in- 
dicate that indeed significant mass growth is likely to be oc- 
curring in the ICL compared to the BCG. Unfortunately the 
semi-analytic models of De Lucia & Blaizot (2007), based on 
the Millennium Simulation and used for comparison in the 
BCG work above, do not tag the stripped stellar content in 
the ICL. More recent work by Puchwein et al. (2010) empha- 
sises just this point. They have carried out high-resolution 
hydrodynamical simulations of clusters, including all known 
physical and feedback processes, to predict the relative build 
up of the ICL and central cluster galaxies. They find persis- 
tently high stellar fractions ~ 40% in the ICL at late times, 
with the vast majority of stars in the main cluster haloes 
belonging to the ICL rather than the BCG - a result which 
underlines the necessity of including both the BCG and ICL 
components when considering the stellar mass budget at the 
centres of clusters. 

The results of Gonzalez et al. (2007) for a sample of 
24 clusters and groups at z ^ 0.13 show an inverse correla- 
tion between the stellar ICL+BCG fraction and the cluster 
velocity dispersion, indicative of the somewhat counter intu- 
itive result that the ICL component grows less efficiently in 
more massive host environments. A non-parametric Spear- 
man rank correlation test applied to our data at the surface 
brightness limit of /ij = 21 mag/arcsec 2 gives a coefficient of 
0.37, with a significance of deviation from zero of 0.47, and 
at fj,j = 20 mag/arcsec 2 we find a Spearman rank coefficient 
of —0.26, with a significance of deviation from zero of 0.62. 
A null result is probably not surprising given our relatively 
small ICL+BCG fractions, although Zibetti et al. (2005) 
also conclude that the ICL fraction remains fairly constant 
(to within ±5%) between different cluster richnesses for their 
SDSS based sample of 683 clusters between 2=0.2 - 0.3. This 
possibly indicates selection differences between the samples, 
however the matter has yet to be settled. 

In our analysis we do not aim to distinguish between the 
ICL and the haloes of the cluster galaxies, but adopt a prag- 
matic ICL definition based simply on an isophotal threshold. 
However, as already mentioned, it is clear in Figures [4] - [9] 
that the extended component we detect closely follows the 
cluster galaxies, especially at higher surface brightness lev- 
els (e.g., fij =19 mag/arcsec 2 in Figures [4]- [9]), rather than 
constituting a diffuse component spread uniformly between 
the galaxies. Our data is not deep enough to detect the pres- 
ence of a faint diffuse (uniformly spaced) ICL component, 
indeed our data would not detect the redshifted diffuse ICL 
component in Virgo at fiv 27 mag/arcsec 2 for example 
(Rudick et al., 2010). As colour gradients in the outer parts 
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of BCGs are shallow, dynamical studies using integral field 
spectroscopy on 8 metre or larger telescopes may be the 
most profitable way forward to separate the BCG halo and 
ICL components (Dolag et al., 20f0). 



6 CONCLUSIONS 

We have detected and measured the fraction of cluster light 
that is in the ICL for 6 galaxies at z ~1 using a simple 
definition of surface brightness threshold. We find that an 
extended component is detectable down to surface bright- 
ness levels of (i j ~22 mag/arcsec 2 measured within a ra- 
dius R500 ■ At this level, the fraction of total measured clus- 
ter light in the ICL for our 6 clusters ranges between 1% 
and 4%, which is smaller than observations at lower redshift 
and similar or slightly below the predicted values at fainter 
isophotes, based on a similar ICL definition from the simu- 
lations of Rll. This indicates that the ICL may have grown 
by a factor of 2 - 4 since z ~1, a scenario which is consis- 
tent with the idea of material being stripped from galaxies 
through mergers and close galaxy encounters. 

In the context of the cosmological mass assembly prob- 
lem of BCGs reported in the literature, the quantity of 
extended light is comparable to the centrally concentrated 
light from the BCGs. Taking into account both components 
is likely to ease the current discrepancy between the ob- 
served and predicted timescales of BCG assembly. 
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APPENDIX A: SURFACE BRIGHTNESS 
PROFILES OF INDIVIDUAL BCGS 
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Figure Al. The ID surface brightness profiles of the BCGs in our sample with best fitting Devaucouleurs profiles. The best fitting 
profiles are convolved with the observed PSF for each fit. 



